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Abstract 

Symplectic and symmetry analysis for studying MHD superfluid flows is devised, a new version 
of the Z. Peradzynski |8] helicity theorem based on differential-geometric and group-theoretical 
methods is derived. Having reanalyzed the Peradzynski helicity theorem within the modern 
symplectic theory of differential-geometric structures on manifolds, a new unified proof and a 
new generalization of this theorem for the case of compressible MHD superfluid flow are proposed. 
As a by-product, a sequence of nontrivial helicity type local and global conservation laws for the 
case of incompressible superfluid flow, playing a crucial role for studying the stability problem 
under suitable boundary conditions, is constructed. 
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1 Introduction 



Long ago it was stated [71 [5] that quantum vortices in superfluid helium can be studied either 
as open hues with their ends terminating on free surfaces of waUs of the container or as closed 
curves. Nowadays the closed vortices are treated as topological objects equivalent to circles. 
The existence of structures such as knotted and linked vertex lines in the turbulent phase is 
almost obvious [12] and forces researchers to develop new mathematical tools for their detailed 
investigation. In this proposed direction it was proved by Z. Peradzyhski [8] a new version of the 
Helicity theorem, based on some differential-geometric methods, applied to the description of the 
collective motion in the in-compressible superfluid. The Peradzyhski helicity theorem describes 
in a unique way, both the superfluid equations and the related helicity invariants, which are, in 
the conservative case, very important for studying the topological structure of vortices. 

Having reanalyzed the Peradzyhski helicity theorem within the modern symplectic theory 
of differential-geometric structures on manifolds, we propose a new unified proof and give a 
magneto-hydrodynamic generalization of this theorem for the case of an incompressible superfluid 
flow. As a by-product, in the conservative case we construct a sequence of nontrivial helicity type 
conservation laws, which play a crucial role for studying the stability problem of superfluid under 
suitable boundary conditions. 

2 Symplectic and symmetry analysis 

We consider a quasi-neutral superfluid contained in a domain M C and interacting with 
a "frozen" magnetic field B : M — > E^, where := (M^, < ., . >) is the standard three- 
dimensional Euclidean vector space with the scalar < .,. > and vector "x" products. The 
magnetic field is considered to be source-less and satisfying the condition B = \7 x A, where 
A : M — > is some magnetic field potential. The corresponding electric field E : M — > E^, 
related with the magnetic potential, satisfies the necessary superconductivity conditions 



where u : M — > T{M) is the superfluid velocity. 

Let dM denote the boundary of the domain M. The following boundary conditions (n, u) \qm = 
and {n,B)\QM = are imposed on the superfluid flow, where n G T*{M) is the vector normal 
to the boundary dM, considered to be almost everywhere smooth. 

Then in adiabatic magneto-hydrodynamics (MHD) quasi-neutral superfluid motion can be 
described, using (12. ip . by the following system of evolution equations: 



E + ux B = 0. 



dE/dt = V X 5, 



(2.1) 




where p : M 



IR+ is the superfluid density, P : M — > E'^ is the internal pressure and 



V 



M 



M is the specific superfiuid entropy. The latter is related to the internal MHD 
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superfluid specific energy function e = e{p, rj) owing to the first thermodynamic law: 



T dr] = de{p, rj) - Pp'^dp, (2.3) 

where T = T{p, rf) is the internal absolute temperature in the superfluid. The system of evolution 
equations ()2.2p conserves the total energy 



H 



^|/i|2 + /5e(/5,r/) + ^|S|2 



d^x, (2.4) 



IM 

called the Hamiltonian, since the dynamical system (j2.2p is a Hamiltonian system on the func- 
tional manifold M := C°°{M; T*{M) xR"^ x E^) with respect to the following [4j Poisson bracket: 



{/> 9} := Im { (M, [|, |].) + P ((|, V|) - (1^ , V|)) 
'i9_M. _ iL^Rw _i_ /R lis. i/i 



+^(V,(gg-S)) + (5,[g,SL) (2.5) 



where we denoted hy p, := pu £ T*{M) the specific momentum of the superfluid motion and by 
[., the canonical Lie bracket of variational gradient vector fields: 

for any smooth functionals f,g £ 'D(M) on the functional space Ai. Moreover, as it was stated 
in [3], the Poisson bracket (|2.5p is, in reality, the canonical Lie~Poisson bracket on the dual space 
to the Lie algebra Q of the semidirect product of vector fields on M and the direct sum of functions, 
densities and differential one-forms on M. Namely, the specific momentum p = pu £ T*{M) is 
dual to vector fields, p is dual to functions, r] is dual to densities and B is dual to the space of 
two-forms on M. Thus, the set of evolution equations (j2.2p can be equivalently re-written as 
follows: 

du/dt = {H, u}, dp/dt = {H, p}, 

(2.7) 

d7]/dt = {H, 7]}, dB/dt = {H, B}. 
The Poisson bracket (j2.5p can be re-written for any f,g£ T>{M) as 

{f,g} = {Df,i^Dg), (2.8) 

with Df := ff^,7^,f^,x4)^ G T*{M) and i9 : T*{M) — > T{M), being the corresponding 



(modulo the Casimir functionals of bracket (j2.5p ) invertible [3] co-symplectic operator, satisfying 
the standard [lOl E] properties 

r = -^9, 5{6w, A i}"^6w) = 0, (2.9) 

where the differential variation complex condition 5^ = is assumed, the differential variation 
vector Sw := {6p,6p,6r],6By £ T*{M.) and the sign "*" denotes the conjugate mapping with 
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respect to the standard bi-linear convolution (., .) of two spaces T*{A4) and T{A4). Note here 
that the second condition of (|2.9p is equivalent [21 [10] to the fact that the Poisson bracket (j2.5p 
satisfies the Jacobi commutation condition. Thus, one can define the closed generalized variational 
differential two-form on A4 

J^^ := {5w,A^-^ 6w), (2.10) 

being a symplectic structure on the functional factor manifold M (modulo the Casimir functionals 
of bracket (f23]l ). 

Denote now a subgroup T>t{M) = {ipt : M — > M} of the diffeomorphism group Diff^M, 
consisting of invertible transformations ipt : M ^ M, generated by MHD superfluid evolution 



equations ()2.2p . This means, by definition, that 

dyDt{x)/dt:=u{^t{x)) (2.11) 
for all X £ M and suitable t G M, for which solutions to p.2p exist and are unique. The symplectic 



structure (I2.10p is invariant with respect to the induced mapping of diffeomorphisms ipt : A4 ^ A4 
on the functional manifold M, that is 

(^^,,..(2) ^ J2) (2.12) 

for suitable t € M. Then the corresponding diffeomorphism subgroup 'Dt{Ai) := {(ft ■ M Ai} 
satisfies the evolution equation 

dipt{w)/dt := Kh{Mw)) (2.13) 

for any w € A4 and the same suitable t G R, where the vector field Kh ■ M — > T{JV[) coincides 
with the system of MHD evolution equations ()2.2I) . This fact easily follows from the standard [2] 
differential-geometric considerations related to equality (12.12p . Really, from (I2.12P one obtains 
that 

= j^^t,.u:^^^ ■■= LKn^^^^ = {iKj + SiK^)J^^ = SiK^co^^^ (2.14) 
for all these suitable t € M, where we denoted by Lk^ the standard Lie derivative with respect 
to the vector field Kh on ^A and used the corresponding Cartan formula = ^Kh^ + SIkh- 
Now, owing to the Hamiltonian equations (|2.7p . the equality ii^-^w^^^ = —5H holds, and since 
5^ = the invariance property (j2.12p is stated. 

As the properties of equations (j2.2p on the manifold M are completely determined by the 
diffeomorphism subgroup 'Dt{M) C Dif f^{M), we will reformulate further the set of equations 
(j2.2p making use of the suitable invariant properties on the manifold M. First, observe that the 
mass conservation law of our superfluid flow is equivalent to the equality 

d 

m 

for any domain Dt CI M moving together with chosen inside particles. It is an easy calculation 
to rewrite (j2.15p in the following equivalent form: 

/ {d/dt + Lu){pd'^x) = Q (2.16) 

JDt 



pd^x = (2.15) 



Dt 



for all domains Df C M and suitable t G M, where as above, we denoted by L„ = i^d + di^ the 
Lie derivative along the vector field u : M — > T[M) on M in the Cartan form. 

As a result of (j2.16p one obtains the following local differential-geometric relationship: 

{d/dt + Lu){pd^x) = f). (2.17) 

Since the evolution of our superfluid is locally adiabatic, the following equality 

{d/dt + Lu)r} = (2.18) 

is obvious, meaning only that dtj/dt = for all suitable t € R. 

Now take the momentum conservation law in the integral Amper-Newton form 

^ [ pud^x+ [ PdSt- [ (j X B) d^x = 0, (2.19) 

JDt JSt=dDt JDt 

where dSt is the corresponding oriented surface measure on the boundary St '■= dDt of a domain 
Dt C M, P : M — > R is the internal pressure and j : M — > is the corresponding induced 
current density in the MHD superfluid under the superconductivity condition. The latter means 
that, owing to neutrality of the superfluid, the induction condition 

VxB + j = (2.20) 

holds. Then from (|2.19|) and (j2.20p one easily obtains the infinitesimal form of the evolution for 
the velocity vector u : M — > T{M): 

{d/dt + Lkjj)u = -p-^VP + p-\V xB)xB (2.21) 



coinciding, evidently, with the first equation of system (j2.2p . 

Consider now at each moment of t E R the subgroup of diffeomorphisms Vr = {ipr : M ^ M} C 
Dif f{M), generated by the following vector field v : M — > T{M) on M: 

diJrix)/dT := viMx)) = p-'B{Mx)), (2.22) 

defined for a suitable evolution parameters r € M. Since the subgroup Vr does not depend 
explicitly on the evolution parameter t G R, its action can be interpreted as re-arranging the 
superfluid particles within any chosen domain Dt C M. Owing now to the commutation property 

[d/dt + Lu,L^]=0, (2.23) 

equivalent to commuting subgroup Dt and Vr for any suitable t, r G R, from the invariance 
condition 

dp/dr = 0, (2.24) 

we can derive that quantities 

7n := L^l (2.25) 
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for all n € Z+ are invariants of the MHD superfluid flow (j2.2p . if the density 7 G A^(M) is also 
an invariant on M. Really, we have 

{d/dt + Luhn = {d/dt + Lu)L^j = K{d/dt + L„)7 = 0, (2.26) 

since, by definition, we have 

{d/dt + Lu)i = 0. (2.27) 

Such a density can be found, observing [4J that the superconductivity conditions E+uxB = 0, 
E = —dA/dt and the last equation of system (j2.2p brings about the invariance condition 

{d/dt + L„)(ia(^) = 0, (2.28) 

where the one-form a^^^ € A^(M) equals 

a^^^ := {A,dx). (2.29) 

Moreover, since the differential operations d/dt + and "d" commute [2], one checks that the 
stronger cohomological condition 

{d/dt + L„)a(^) = (2.30) 

holds on M, if the time-dependent gauge mapping A — > A+Vip, where dip/dt+Luip+{u, A) = 0, 
is applied to the magnetic potential A : M — > E^. Now from conditions (j2.28p and (j2.3Up one 
easily derives that the density 

^ := a(i) A (ia(i) (2.31) 

satisfies equation ()2.27p . Thus, it generates, in view of formula (j2.25p . new conserved quantities, 
which can be equivalently rewritten as 

7„ := pL:{p-\B,A)) = pL:{v,A) (2.32) 

for all n G Z+. Thereby, the following functionals on the functional manifold A4 

Hn.= [ lnd^x= [ pL^{p-\B,A))d^x (2.33) 
Jm Jm 

for all n G are invariants of our MHD superfluid dynamical system (|2.2p . In particular, at 
n = we obtain the well-known [4J magnetic helicity invariant 

Hq= [ {A,V X A) d^x, (2.34) 

which exists independently of boundary conditions, imposed on the MHD superfluid flow equa- 
tions ([2:2|) . 

The result obtained above can be formulated as the following theorem. 

Theorem 2.1 The functionals \2.33\) . where the Lie derivative Ly is taken along the magnetic 
vector field v = p^^B, are global invariants of the system of compressible MHD superfluid and 
superconductive equations 112. 
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Below we proceed to symmetry analysis of the incompressible superfluid dynamical system 
and construct the related local and global new helicity invariants. The case of superfluid hydrody- 
namical flows [9j is of great interest for many applications owing to the very nontrivial dynamical 
properties of so-called vorticity structures, featuring the motion. 

3 The incompressible superfluid: symmetry analysis and conser- 
vation laws 

Concerning the helicity theorem result of [8j, where the kinematic helicity invariant 

Hq := / (u, V X u) (fx (3.1) 

JM 

was derived, making use of differential-geometric tools in Minkowski space in the case of incom- 
pressible superfluid at the absent magnetic field i? = 0, we will show below its general dynamical 
symmetry nature. The governing equations look as follows: 

du/dt = -{u,V)u + p-^VP, dp/dt + {u,Vp) =0, (V,n) = 0, (3.2) 

where the density conservation properties 

{d/dt + Lu)p = 0, {d/dt + Lu)d^x = (3.3) 

hold for all suitable i G M. Define now the vorticity vector ^ := V x u and find from (13. 2p that it 
satisfies the vorticity flow equation 

a^/at = V X (n X ^). (3.4) 

Really, the first equation of (|3.2|) can be rewritten as 

du/dt = n X (V X ti) - p'^VP - ^V|ttp. (3.5) 

Then, applying the operation "V x • " to (|3.5p . one easily obtains the vorticity equation (13. 4p . 
Moreover, equation (j3.4p can be rewritten in the equivalent form 

di/dt+{u,V)i = {i,V)u, (3.6) 

which allows a new dynamical symmetry interpretation. 
Put, by definition, 

dx/dr = v{x,t) := p'^i, (3.7) 

defining for all r € M the diffeomorphism subgroup T>r C Diff M of the manifold M. It is easy 
to check that this subgroup commutes with the previous defined subgroup Vf C Diff M, since 
the following condition 

{d/dt + Lu)v = L^u (3.8) 
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holds for all t,T € M, exactly coinciding with relationship (j3.6p . The condition (jc{.8|) means that 
the commutation property 

[9/9t + L„,L^] =0, (3.9) 

similar to (j2.23p . holds. 

Now we can make use of the invariants generation technique, described above in the case of 
the superfluid equations (j2.2p . For this we need to construct a source density invariant 7 € (M) 
of equations (|2.26p and construct successively a hierarchy of additional invariants as 

In ■■= L^.l (3.10) 

for all n G 

Put, by definition, G K^{M) as the one-form 

:= (u,dx) (3.11) 

and find that 

{d/dt + L„)/3W = -p-^dP + ^d|n|2 = dip'^P + ^|u|2). (3.12) 
The differential two-form d/J^^-* G A^(M) satisfies the condition 

{d/dt + Lu)dp^^^ = d^{p-^P + ^|up) = (3.13) 

owing to the identity d^ = 0. Then the differential density three-form 7 := [3^^^ A d/?^"'^-* G h?{M) 
satisfies, owing to (|3.12p and (|3.13p . the condition 

{d/dt + L„)7 = {d/dt + Lu){f3^^^ A d/3(i)) 

(3.14) 

= d {p-^P + i|n|2) A = d{{p-^P + i|n|2) d/3W) . 

By integration of (j3.14p over the whole manifold M we obtain, based on the Stokes theorem, the 
expression 



I /3(y A = i !^,{u X (V X u)) d^x = i f^,{n x ^ d'x 

(3.15) 



if the boundary conditions {u,n) = and S,\dM = are imposed on the superfluid vorticity flow. 
Really, the surface measure {du, Adx) on the boundary dAI can be equivalently represented as 

{du, Adx) = {{dx, V)u, Adx) = (V x u, dS) = dS), (3.16) 

where dS is the standard oriented Euclidian surface measure on dM. Since the vorticity vector 
CldM = 0, the result (|3.15p follows automatically. 
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Assume now that the vorticity vector ^ = V x u satisfies the additional constraints L'^^\dM = 
for n € Z_|_. Then we obtain from (j3.16p and (|2.23p that 

ft hi ^"7 = i Im ^"7 d'^ = J^, md/dt + L„)7 
= f,,Ly{p-'P+l\u\') dpW = J^^dL^ i[^,|2) dpW) 

= Jj^dL- + {dn,Adx)) = f^.dL^ ((p-ip + i |n|2) {^,dS)) 

= !,^,L-, (e,(i5)) = U,Y.l=,C':,{L%^-^ 1|^12) dS)) = 0, 

(3.17) 

bringing about the generahzed hehcity invariants 

Hn:= [ pL'^,{uxC)d'^x (3.18) 
Jm 

for all n € Z^. Notice here that all of the constraints imposed above on the vorticity vector 
^ = V X n will be automatically satisfied, if the condition supp ^ fl dM = holds. The result 
obtained can be formulated as the following theorem. 

Theorem 3.1 Assume that an incompressible superfluid, governed by the set of equations /i3. 2\) 
in a domain M C E^, possesses the vorticity vector ^ = V x n, which satisfies the boundary 
constraints L'^_i^S,\dM for all n E Then all functionals h3.1^) will be the generalized helicity 
invariants of \3.2(l . 

The results obtained above allow some interesting modifications. To present them in detail, 
observe that equality (j3.12p can be rewritten as 

{d/dt + L„)/3(^) -dh = (d/dt + L„)/3(^) = 0, (3.19) 

where, by definition, 

/i := i|up, := - Vv9,dx), (3.20) 
and the scalar function : M — > M is chosen in such a way that 

{d/dt + Lu)ip = Vh. (3.21) 

Then, obviously, there holds the additional equality 

{d/dt + Lu)d/3^^'> = 0, (3.22) 

following from the commutation property [d,d/dt + Ly] = 0. Then we obtain that the density 
p := A dp'-^^ G A^{M) satisfies the condition 

{d/dt + L„)/i = 0, (3.23) 

for all t G M. The similar equality holds for densities jln ■= L^p & A^(M), n € Z+i 

{d/dt + Lu)fln = 0, (3.24) 
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owing to the commutation property (j2.23p . Thereby, the fohowing functionals on the correspond- 
ing functional manifold A4 are invariants of the superfiuid flow (j2.26p : 

Ttn:= [ fin= [ pL"" iJiu-V^),0 <fx (3.25) 

Jm JDt 

for all n € Z-|_ and an arbitrary domain Dt C M, independent of boundary conditions, imposed on 
the vorticity vector ,^ = V x li on dM. Notice here that only invariants ()3.25p strongly depend on 
the function ip : M — > M, implicitly depending on the velocity vector u E T{M). Mention here 
only that the practical importance of the constructed invariants (|3.25|) is not still clear enough. 

4 Conclusions 

The symplectic and symmetry analysis of compressible MHD superfiuid, as shown above, appeared 
to be effective for constructing the related helicity type conservation laws, important for practical 
applications. In particular, these conservative quantities play a decisive role [HH], when studying 
the stability of MHD superfiuid flows under special boundary conditions. 

Here we also need to notice that the differential-geometric reformulation of MHD equations 
(j2.2p suggested in 0] is incorrect. Namely, the equality {d / dt+ L^j^) {p^^ B , dx) = is not equivalent 
to the magnetic field equation dB/dt — V x (V x i?) = that one can check by easy calculations. 
Nonetheless, the commutator relation [d/dt + L„,Lp-i^] = devised there and all Casimir 
invariants found in article 0] are true. But some problems related tothe construction of non- 
Casimir type MHD superfiuid flows using their Hamiltonian structure remain, in general, open 
and wait still to be treated in detail. Some of the results in this direction can be eventually 
obtained making use of group-theoretical and topological tools developed in [THl lllj . where 
the importance of the basic group of diffeomorphisms Diff{M) of a manifold M C and its 
differential-geometric characteristics were stated. 
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